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Abstract 

We propose a method to determine the locally preferred structure of model liquids. This latter 
is obtained numerically as the global minimum of the effective energy surface of clusters formed 
by small numbers of particles embedded in a liquid-like environment. The effective energy is the 
sum of the intra-cluster interaction potential and of an external field that describes the influence 
of the embedding bulk liquid at a mean-field level. Doing so we minimize the surface effects 
present in isolated clusters without introducing the full blown geometrical frustration present in 
bulk condensed phases. We find that the locally preferred structure of the Lennard-Jones liquid is 
an icosahedron, and that the liquid-like environment only slightly reduces the relative stability of 
the icosahedral cluster. The influence of the boundary conditions on the nature of the ground-state 
configuration of Lennard-Jones clusters is also discussed. 
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I. INTRODUCTION 



A liquid is said to be supercooled when it is possible to cool it below its melting temper- 
ature T m without crystallizing. The supercooled liquid phase is metastable with respect to 
the underlying crystal, and it is characterized by a dramatic increase of the viscosity and 
the relaxation times upon lowering the temperature, an increase that eventually leads to 
glass formation. These dramatic changes in the dynamic properties are not accompanied 
by strong signatures in the structural quantities, such as the static structure factor. Yet, 
it has been suggested that both supercooling and glass formation were deeply connected 
to the structure of the liquid, more precisely to a competition between extension of a local 
liquid order, different than that of the crystal, and global constraints associated with tiling 
of the entire space 

HQS. 

This competition has been termed geometric (or topological) 

frustration |jl Q • 

Some fifty years ago, Frank put forward the following argument to explain supercooling 
of liquids |4|. If one considers atomic liquids in which atoms interact via spherically sym- 
metric potentials like the Lennard- Jones potential, the local arrangement of the atoms that 
is preferred is not the clusters associated with crystalline order (face-centered cubic and 
hexagonal close-packed lattices), but a polytetrahedral packing, the icosahedron formed by 
13 atoms. The energy of such an icosahedral arrangement interacting via the Lennard- Jones 
potential is indeed 8.4% lower than the close-packed crystalline clusters. Local icosahedral 
order should then be prevalent in liquids but, because of the 5-fold rotational symmetry 
of the icosahedron, it cannot tile the entire space and form a crystal: this is a manifesta- 
tion of geometrical frustration. Crystallization then requires a rearrangement of the local 
structure of the liquid, which leads to a strong first-order freezing transition and allows 
supercooling of the liquid. Since then, there has been a large body of experimental and 
simulation work confirming the prevalence of local icosahedral, or more gen erall y p olvtetra- 
hedral, order in atomic liquids and metallic glasses [3, 0, S fl 0, 0, Qllll Q Q- 

The tendency to form icosahedral order has been shown to increase as the temperature is 



lowered 



dency tc 
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Frank's argument has, however, one shortcoming 16]: the 13-atom cluster considered is 
isolated, so that most of the energetics is related to the surface, a situation that of course 
does not occur in bulk liquids. How can one remedy this problem? An a priori easy way 
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would be to study directly the local arrangement of the atoms in a bulk liquid. In such a 
case, however, it is found that the proportion of icosahedra among all the 13-atom groups 
is very small, typically a few percent [10]. This is indeed to be expected since geometrical 
frustration is present, which opposes the growth of icosahedral order and distorts the local 
polytetrahedral arrangements P, How to disentangle then the determination of the 
locally preferred structure from frustration effects? The method we propose in this work 
is to consider the influence of the bulk liquid on a given 13-atom cluster at some mean- 
field level, so that surface effects can be reduced and made more realistic for describing a 
condensed phase, whereas geometrical frustration is strongly inhibited. The main advantage 
of this method is that it can be extended to study the locally preferred structure of molecular 
liquids, for which a priori topological arguments do not easily provide the symmetries of all 
possible local arrangements of the molecules, nor the nature of geometrical frustration. 

More specifically, in this article we consider the ground-state of a cluster of 13 atoms in- 
teracting via a Lennard- Jones potential (smaller and larger clusters are also considered); the 
atoms are placed in a cavity and are subject to an external field that mimics the interaction 
with the rest of the liquid. The structure of the outside liquid only enters the calculation 
via the bulk pair distribution function (known from previous simulation studies). By means 



of an optimization algorithm 



18, 



19( we find that the ground-state of the cluster, i.e., the 



global minimum of the (effective) energy surface formed by the intra-cluster interactions 
and the external field, is of icosahedral symmetry, therefore generalizing Frank's result. For 
sake of comparison we consider in addition other boundary conditions for the cluster. These 
conditions describe different types of environments: free boundary conditions (isolated clus- 
ter), periodic boundary conditions (periodic tiling of space), and icosahedral-like boundary 
conditions (hypothetical non-frustrated system). 



II. METHOD AND CHOICE OF THE BOUNDARY CONDITIONS 

We consider a system of atoms interacting via a pair-wise additive spherically symmetric 
potential v(r), where r is the center-to-center distance. A number N of atoms are placed in 
a spherical cavity C of radius Rc, that we envisage as surrounded by bulk liquid made of the 
same atoms, and characterized by the temperature T and the density p. As explained above, 
we do not want to fully account for the liquid structure because geometrical frustration would 
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obscure the nature of the local order. We rather resort to a mean-field type of description in 
which the liquid outside the cavity is considered as a continuum, characterized by a (known) 
pair distribution function g(r;T,p) that is not affected by the fact that a cavity has been 
carved out. The potential energy acting on a given atom at position r inside the cavity due 
to the outside liquid is thus described as 

W{v-R c ) = p - [ d 3 v> sdr-r'IMIr-r'l), (1) 

where the integral is over all positions outside the cavity. Taking the center of the cavity as 
origin and transforming to spherical coordinates, one finds after standard manipulations 

{[•Rc+r i-oo \ 

I dx x 2 g(x)v (x) [1 — u(x; r, Rc)} + 2 / dx x 2 g(x)v (x) > , (2) 
JRc-r JRc+r ) 

where 

u(x;r,R c )=( Rh ~ r2 ~ x2 ). (3) 
\ 2rx J 

Note that u{x = Rc — r) — —u(x = Rc + r) = 1. The detail of the calculation is given in 
the Appendix „ 

In view of the implementation of the optimization algorithm we need an explicit 

expression for the first derivatives of the external potential, both with respect to Rc and r; 
we obtain 

Rc+r 
Rc-r 

Rc+r f R% + r 2 — X 2 

dx x 2 g(x)v(x) ( 

R C -r 

The second derivatives, needed for the calculation of the Hessian matrix and the study of 
the tuition states 00,0^ be ob^ed in a si^wa, 

The total potential energy for the iV atoms of the embedded cluster is the sum of the 
atom-atom interaction potentials inside the cavity and of the external potential, 

N N 

"({'•,}, v :/,v ) = E <\ r i- r i\)+J2 w ( r fi R ^ ( 6 ) 

i>j=l j=l 

where W(r; Rc) is given by Eqs. (0) and (jSJ). Finding the ground-state configuration for 
the iV-atom cluster embedded in a liquid-like environment (at a given T and p) amounts 
to determining the global minimum of U with respect to variations of the positions of the 
iV atoms. (Note that T and p only enter trough the external potential W(r;R c ), both 



d 
dRc~ 

and 



W(r;R c ) =-^p(—^\ \° dx xg(x)v(x), (4) 

r \ r J J Rc-r 



d_ 

dr 



W{r; R c ) = up I dx x 2 g{x)v{x) c - ' - - ) . (5) 
Rc J Rn -r V 2xr 2 



explicitly (see Eq. (J2J)) and through the state-dependence of the pair distribution function 
g(r).) It is worth stressing that, contrary to the case of an isolated cluster for which the 
radius of the cavity Rc is merely fixed to avoid evaporation of the atoms, Rc becomes 
a relevant variable in a liquid-like environment: to preserve a realistic description of the 
liquid, Rc should adjust to global contractions or expansions of the iV-atom cluster, which 
should then be taken into account in the minimization procedure. This point will be further 
discussed below. 

To obtain the lowest energy minimum of U, we have used a slightly modified version of 
the basin-hopping algorithm introduced by Wales and co-workers jl8L Iioj]. The algorithm 
consists of a constant-temperature Monte-Carlo simulation performed with an acceptance 
criterion based not upon the energy of the proposed new configuration, but upon the energy 
of the closest minimum of the potential energy surface, obtained by a local minimization 
starting from that configuration. This algorithm turns out to be a very efficient method for 
exploring directly the minima of the potential energy surface, and it allows one to locate the 
ground-state with relatively little effort. 

Finally we have also considered other boundary conditions for the iV-atom cluster: 
i) Free boundary conditions (isolated cluster).- This is the standard case studied in the 
literature. It simply corresponds to the above situation in which the external potential is 
set to zero: 



ii) Periodic boundary conditions (periodic replication of the local cluster).- Each atom of 
the cluster now interacts also with the images of the other atoms of the cluster 



where Sj is the position of the image of atom j selected through the minimum image criterion, 
i.e., among all possible images of atom j, only the closest is selected. We have used a cubic 
elementary cell. 

Hi) Icosahedral-like boundary conditions (hypothetical non-frustrated system).- As men- 
tioned in the Introduction, icosahedral order cannot be extended to the entire space. One can 
however introduce icosahedral-like boundary conditions by embedding the 13-atom cluster 
in the center of a large 147-atom cluster with icosahedral (J/j) symmetry in which the central 
atom and its first layer have been removed. (Recall that atomic clusters are characterized 



W(r;R c ) = 0. 



(7) 




(8) 
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by "magic numbers" of atoms for which the global minimum has polytetrahedral symme- 



try 



20 



211 ]: in particular, icosahedral symmetry is obtained for 13, 55, and 147 atoms.) The 



external potential is now written as 

W(r)= 1 - ^(l r - r *D> (9) 

he Qico qico 

rii tZ '-'2 ? 3 

where S^ 00 and are the second and third shells of the 147-atom icosahedron. For com- 
pleteness, we have also considered the 55-atom Mackay icosahedron, as well as the correction 
due to embedding the 55-atom or 147-atom icosahedron in bulk liquid with the resulting 
external potential treated at a mean-field level (see above). 



III. GLOBAL MINIMUM OF THE LENNARD-JONES CLUSTERS 



We specialize our investigation to the case of the Lennard- Jones pair potential 



vlj{t) = 4e 



0-N !- 

r . 



<7\" 

r 



(10) 



where e and 2 1 / 6 cr are the well-depth and the separation at the minimum of the potential, 
respectively. In what follows we set a = e = 1. 

In order to evaluate the liquid-like external potential acting on the iV-atom cluster, 
Eqs. (J2J and ©, and its derivatives, Eqs. (J1J and ©, one needs a model for the pair 
distribution function g(r). We use the 7-parameter parametrization of Verlet's Molecu- 



lar Dynamics simulation data on the Lennard- Jones liquid 
Mansoori 



221 ] proposed by Matteoli and 



g(y) = l + y~ m [g(d) 
for m > 1, y > 1, and 



1-A] + 



(y 



A) 



exp [-a(y - 1)] cos[/% - 1)], 



'ID 



g(y)=g(d)exp[-6(y-l)% (12) 

for y < 1. Here y = r/d is the dimensionless intermolecular distance where d = 2 1 / 6 cr, 
and h, m, A, a, /3, 6, g(d) are adjustable parameters. The terms y~ m and exp [— 8(y — l) 2 ] 
describe the decay of the first peak, while the term exp [—a(y — 1)] cos[/3(y — 1)] provides the 



23j. We have taken the values h = 1.065, 



damped oscillations observed at larger distances 
m = 13.42, g(d) = 2.830, A = 0.9310, a = 1.579, (3 = 6.886, 6 = 135.9 that allow to 
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reproduce the pair distribution function for the liquid at p = 0.880 and T = 1.095 (in usual 
reduced Lennard- Jones units) 



22j,|2i 

The resulting external potential W(r; Rc) is shown in Fig. Q for several values of the 
cavity radius Rc- The shape of the r-dependence changes with Rc so that no rescaling of 
the curves is possible. For Rc < 1.2 a, the potential has a minimum at r = because the 
central atom sits at the minimum - or very close to it - of the pair interactions due to liquid 
atoms at the boundary of the cavity; this is no longer true for larger cavity radii, and W 
decreases monotonically with r, the most favorable position inside the cavity being at its 
edge where the attractive interaction due to the nearby liquid particles is the strongest. By 
construction, when R c — > 0, W(r = 0, Rc) becomes equal to the total potential energy of 
the Lennard- Jones liquid at the considered state point, E ~ —5.7 e, whereas when Rc — > oo, 
W(r — Rc, Rc) is equal to half this energy. 

As we discussed it above, the radius of the cavity Rc must be adjusted to global contrac- 
tions or expansions of the cluster. A reasonable way to implement this is to take at each 
minimization step, i. e., for each configuration of the N atoms, 

Rc = r max + /i a, (13) 

where r max is the distance of the outermost atom from the center of the cavity, and /z is a 
constant chosen to account for the fact that repulsive interactions between atoms make very 
unlikely the presence of "bulk liquid" atoms when their centers are too close to those of the 
cavity atoms; we have taken /i = 0.5, but we have checked that the results are independent 
of the actual value, in runs with different values of /z between 0.1 and 1.0. 

A typical minimization run for a 13-atom cluster in the presence of a mean-field liquid-like 
environment is shown in Fig. |21 where we have plotted the evolution of the energy and its 
intra-cluster and external- field contributions (bottom), together with the evolution of the 
cavity radius Rc (top). During the optimization run that starts from a random configuration 
of 13 atoms, the structure of the cluster becomes more and more compact, and its energy 
decreases. The final optimized configuration is found to have icosahedral symmetry. Both 
the final cluster radius, r max ~ 1.08 a, and the final intra-cluster energy, Ui ntr a — —44.327 e 
are identical to those found for the ground-state of the isolated 13-atom cluster j^J. The 
total energy of the icosahedral cluster in the presence of a liquid-like environment is however 
much lower (by almost a factor two) because of the external field that compensates for the 



deficiency of nearest neighbors of the 12 surface atoms. We have repeated the procedure for 
different starting random configurations and we have always obtained the icosahedron with 
r max ~ 1.08 a as the global minimum. The same is true for a whole range of liquid density 
p and temperature T (0.65 < p < 0.88, 0.6 <T< 3.6). 

In addition to carrying out a global optimization procedure, one may also make a calcula- 
tion in the spirit of Frank's pioneering work |4] . We have considered two potential candidates 
for the ground-state configuration of the 13-atom cluster, namely the icosahedron and the 
cuboctahedral cluster with Oh symmetry that is associated with the face-centered-cubic 
close-packed lattice and we have compared their energies in a liquid-like environment. The 
energies of the two, /^-symmetric and O/j-symmetric, clusters are shown in Fig. |H]as a func- 
tion of the cavity radius Rc (or, equivalently, as a function of r max , a unique distance being 
enough to fully determine the whole cluster once the symmetry, J/j or Oh, is chosen.). The 
intra-cluster contribution to the energy has a minimum for r max ~ 1.08 o for the icosahedron, 
and r max ~ 1.10 a for the Oh cluster, whereas the external-field contribution monotonously 
increases with Rc in both cases. One then finds that the icosahedral cluster has a minimum 
total energy for r max ~ 1.08 a and that this energy is 4.8% lower than the lowest energy 
found for the Oh cluster when r max ~ 1.10 a. When compared to Frank's result (8.4%), the 
relative energy difference between the two types of clusters is thus not drastically modified: 
a mean-field liquid-like environment only slightly reduces the relative stability of the icosa- 
hedral order. Finally, as a mere check of our global optimization procedure, we have verified 
that the icosahedral ground-state found here is identical to that discussed above. 

The influence of the boundary conditions on the ground-state of a 13-atom cluster can be 
investigated by using again the optimization algorithm (always starting with a random initial 
configuration) with the appropriate conditions, free, periodic, and icosahedral-like, described 
in the previous Section. As already well-known, the ground-state of the isolated cluster 
is an icosahedron and, as anticipated, that of the cluster in the presence of icosahedral- 
like boundary conditions is also an icosahedron. The ground-state energies are shown in 
Table HJ and can be compared to that of the icosahedral cluster in a mean-field liquid-like 
environment. Not surprisingly, this latter is much lower than that of the isolated cluster (see 
above), but it is higher than that of icosahedra embedded in larger icosahedral structures. 

We note on passing that, in the case of an icosahedral-like environment, the change in 
the structure of the 13-atom cluster from random to icosahedral during the optimization 
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run starts from the outside and propagates inward. This is reminiscent of what has been 



observed in the simulation of gold nanoclusters [2J|. There, it has been found that, just after 
freezing, ordered nanosurfaces with five-fold rotational symmetry are formed, while interior 
atoms remain in a disordered state. On lowering the temperature, the crystallization of 
the interior atoms proceeds from the surface toward the core region, eventually producing 



an icosahedral structure [2J| . This is at variance with the classical picture of homogeneous 
nucleation and rather represents a surface-induced (heterogeneous) crystallization. 

Periodic boundary conditions lead to a quite different picture. Since icosahedra cannot 
tile space by periodic replication, such conditions should favor the symmetries that allow 
a complete filling of space with true long-range order. It is indeed what we have found: 
the global minimum is then a cuboctahedral cluster with Oh symmetry, that leads to a 
face-centered-cubic close-packed lattice when periodically replicated (see Fig. HJ); the cor- 
responding ground-state energy is given in Table HI and it is found lower than that of an 
isolated icosahedron — because of the lack of neighbors already mentioned for this latter 
case — but higher than icosahedra in either liquid-like or icosahedral- like environment. 

Finally, we have considered the effect of varying the number of atoms present in the 
cluster. The main motivation for this study is to check that without a-priori knowledge of 
the preferred local structure in the presence of a mean-field liquid-like environment, hence 
of the number of atoms involved in this structure, the global minimization method will help 
select the proper preferred configuration. This will be important when considering molecular 
liquids. We have thus studied the global minima of iV-atom clusters with N ranging from 
2 to 23. As shown in Fig. El the energy per-atom, the only relevant quantity for comparing 
local structures in a liquid-like environment, is lowest for the N = 13 icosahedral ground- 
state. For small N, the mean-field description of the liquid environment is probably too 
crude to give sensible results, but it is nonetheless significant that for a large range of N, 
the icosahedral cluster is properly selected as the locally preferred structure of the Lennard- 
Jones liquid. 

IV. CONCLUSION 

Local icosahedral order has been found both in bulk condensed phases and in clusters 
formed by spherical particles. The connection between bulk and cluster studies is however 
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obscured by two facts: first, geometrical frustration strongly hinders the spatial extension 
of local icosahedral order and distorts the local icosahedra in bulk conditions; second, the 
energetics of isolated clusters is partly determined by surface effects that are of course absent 
in the bulk. In this work we have tried to bypass those problems. We have proposed to 
determine the locally preferred structure of a liquid by finding the ground-state configuration 
of iV-particle clusters embedded in a liquid-like continuum, characterized by the proper 
density and pair distribution function of the bulk liquid at the chosen thermodynamic state 
point. This mean-field-like procedure minimizes the surface effects without introducing full 
blown geometrical frustration. 

In terms of potential energy surface, we have therefore introduced an effective energy 
surface that contains the usual intra-cluster potential energy contribution plus an external 
field that accounts for the interaction with the outside liquid at a mean-field level. By 
a global optimization algorithm we have then located the lowest-energy minimum of the 
effective energy surface. For Lennard- Jones pair interactions, we have found that the locally 
preferred structure is indeed an icosahedron, the effect of the liquid-like environment being 
to only slightly reduce the relative stability of the icosahedral structure when compared 
to Frank's calculation for an isolated cluster j^. We have also shown the importance of 
the boundary conditions used for the cluster: whereas icosahedral- like boundary conditions 
stabilize even more the local icosahedral cluster, periodic boundary conditions make the 
cuboctahedral cluster more stable, a consequence of geometrical frustration that prevents 
tiling of space by icosahedra, and therefore favors long-range order associated with face- 
centered cubic or hexagonal close-packed lattices. 

The present findings for the Lennard- Jones liquid suggest that the proposed method could 
be efficient as well for determining the locally preferred structure of molecular liquids, in 
cases where both translational and rotational degrees of freedom are involved, and a-priori 
knowledge about the putative local order is scarce. 
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APPENDIX: MEAN-FIELD LIQUID EXTERNAL POTENTIAL 



Here we give some details about the calculation of the mean-field liquid external potential 
W(r; Re) discussed in Sect. HT1 

From Eq. (|1J) the external potential felt by an atom at r, where the origin of the coordi- 
nates is chosen as the center of the cavity, can be written as 

W(r,Re) = £f d 3 xg{x)v(x), (A.l) 

1 J\x+r\>R c 

where we have changed the integration variable from r' to x = r - r'. We now rotate the 
reference system such that r/r = z and we translate it so that the new origin is at P, the 
position of the center of the atom under consideration. The new geometry of the problem 
is shown in Fig. El One can easily convince oneself that 

/*7T /*27T /*00 

W(r;R c ) = - d<f> sin0 / d6 dx x 2 g(x)v(x). (A.2) 
2 Jo Jo J sty) 

Given a point inside the cavity, s is its distance from the surface in the direction x. It is 
simple to check that s is a function of (f> only and is solution of the quadratic equation 

s 2 + (2r cos <j))s + (r 2 - R c ) 2 = 0; (A.3) 

the correct solution is 



= — r cos + y r 2 cos 2 <p + (R 2 C — r 2 ). (A. 4) 

Introducing the variable u = cos we obtain 

/l roc 
du I dx x 2 g(x)v (x), (A.5) 
-1 Js(u) 

with 



s(u) = -ru + yr 2 {u 2 - 1) + R 2 C . (A.6) 

s{u) is such that s(— 1) = Rc + r, and s(l) = Rc — r. We now change variable from u to s, 
and integrate by parts so that 

" Rc ~ r , du r 



W(r; Rc) = np ds — / dx x 2 g(x)v(x) 

Jrh+t ds J ' 



up 



Rc+r us J P 

poo 

u(s) I dx x 2 g(x)v(x) 
11 



p 



Rc-r j-Rc+r 

- I ds s 2 u(s)g(s)v(s) } , (A. 7) 

R C +r JRc-r 



where 



and u(s = Rc — r) — —u(s = Rc + r) = 1. Substituting Eq. (|A.8J) in Eq. (|A.7J) and 
rearranging we obtain Eqs. (j2J) and (jHJ). 
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Boundary Conditions 


Energy [e] 


Point Group 


Free 


-44.327 


h 


Periodic 


-52.745 


o h 


Liquid-like 


-88.490 


h 


Icos. (1 lay.) 


-92.980 


h 


Icos. (2 lays) 


-100.142 


h 


Icos. (1 lay.) & Liquid 


-101.395 


h 


Icos. (2 lays) Sz Liquid 


-102.208 


h 



TABLE I: Ground-state energy and symmetry of the 13-atom cluster with various boundary condi- 
tions. The icosahedral-like conditions correspond to the 55-atom (1 layer) and 147-atom (2 layers) 
icosahedral clusters, and the same structures embedded in a mean-field liquid-like environment. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

r[o] 



FIG. 1: Mean- field external potential as a function of the distance r from the center of the cavity ; 
different curves are for different values of the radius of the cavity Rc- The dotted line is the value 
of the potential at the surface of the cavity. 
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1 1 ' 1 ' 1 ' 1 ' 1 

50 100 150 200 

Minimization Steps 



FIG. 2: Evolution of the cavity radius (Top) and of the cluster energy, together with its intra-cluster 
and external- field contributions (Bottom), during a minimization run starting from a random 13- 
atom Lennard- Jones cluster with mean- field liquid- like environment. 



15 




1.2 1.4 1.6 1.8 1.2 1.4 1.6 1.8 

FIG. 3: Energy of a 13-atom cluster with mean-field liquid-like environment as a function of 
Rc = r max + l~i cr (diamonds): (a) icosahedral (Ih) cluster, (b) f-c-c {Oh) cluster. We show 
separately the contributions of the intra-cluster interaction energy (circles) and that of the external 
field (squares). 




( a ) I h symmetry ( b ) O h symmetry 



FIG. 4: Structure of the icosahedral (a) and cuboctahedral (b) 13-atom clusters. 
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